Electronic transport in Si nanowires: Role of bulk and surface disorder 
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We calculate the resistance and mean free path in long metallic and semiconducting silicon 
nanowires (SiNWs) using two different numerical approaches: A real space Kubo method and a 
recursive Green's function method. We compare the two approaches and find that they are com- 
plementary: depending on the situation a preferable method can be identified. Several numerical 
results are presented to illustrate the relative merits of the two methods. Our calculations of relaxed 
atomic structures and their conductance properties are based on density functional theory without 
introducing adjustable parameters. Two specific models of disorder are considered: Un-passivated, 
surface reconstructed SiNWs are perturbed by random on-site (Anderson) disorder whereas defects 
in hydrogen passivated wires are introduced by randomly removed H atoms. The un-passivated 
wires are very sensitive to disorder in the surface whereas bulk disorder has almost no influence. 
For the passivated wires, the scattering by the hydrogen vacancies is strongly energy dependent and 
for relatively long SiNWs {L > 200 nm) the resistance changes from the Ohmic to the localization 
regime within a 0.1 eV shift of the Fermi energy. This high sensitivity might be used for sensor 
applications. 

PACS numbers: 73.63.-b, 72.15.Lh, 72.10.Fk 



I. INTRODUCTION 

Semiconducting nanowires are very promising building 
block for future nanoelectronic and nanophotonic appli- 
cations as witnessed by several recently demonstrated 
deviceaiSiiii^. Silicon nanowires (SiNWs) are espe- 
cially attractive candidates due to their compatibility 
with conventional Si-technology and due to the accu- 
rate control of diameter and electronic properties dur- 
ing synthesis^. Furthermore, in recent years SiNWs have 
been applied as label-free real-time chemical and biologi- 
cal sensors with very high sensitivity and, e.g., capability 
of single virus detection^. 

Thin SiNWs with diameters below 5 nm have been syn- 
thesized by several groups. Recently, Ma et al^ obtained 
very thin wires grown in the (110) and (112) directions 
with diameters down to 1.3 nm, and Holmes et ali^ have 
previously reported 4-5 nm (100) and (110) SiNWs. Wu 
et al7 recently demonstrated that the growth directions 
depend on the diameter, which can be controlled by the 
size of a catalytic nanoparticleii. 

Concerning theoretical modelling of the structural 
properties, Rurali and Lorente-^ recently showed, us- 
ing ab initio calculations, that thin (100) un-passivated 
SiNWs could be either metallic or semi-metallic depend- 
ing on the surface reconstruction. In another recent work, 
Singh et ali^ theoretically studied pristine (110) SiNWs 
and found that these were indirect band gap semicon- 
ductors. Very recently, Fernandez-Serra et ahM used 
ah initio calculations to study the surface segregation of 
dopants in both passivated and un-passivated SiNWs. 

The large surface to bulk ratio in the thin wires imply 
that surface effects such as defects, vacancies or adatoms 
will have a large influence on the transport properties. 



This was experimentally demonstrated by Cui et alX 
showing increased mobilities after passivation of surface 
defects. Also, electron-phonon scattering might be sup- 
pressed in thin wires. Indeed, recent experiments by Lu 
et alM> indicated ballistic transport in undoped Si/Ge 
core-shell wires at room temperature with an estimated 
phonon scattering mean free path (MFP) Iph > 500 nm. 
This might imply that even at room temperature de- 
fects could be the most important scattering source, and 
a thorough understanding of the scattering processes is 
thus required. 

A number of transport calculations have been reported 
for wires with various diameters. Das et al.^- used the 
Boltzmann equation to calculate the carrier mobility in 
relative thick (d = 10 — 90 nm) GaAs wires focusing on 
the diameter dependence. Sundaram et al'M- also used 
the Boltzmann equation to study surface effects on the 
transport in large diameter wires. Zheng et ali^ applied 
a tight-binding model of a hydrogen passivated wire and 
studied the effect of wire thickness on the band gap, ef- 
fective masses and transmission. 

Real SiNWs with lengths up to the scale of microm- 
eters consist of millions of atoms and are likely to have 
many randomly placed defects. To our knowledge, no 
theoretical works concerning SiNWs, based on ab initio 
methods and including many scattering events, have yet 
been published. 

A calculation of the conductance of a SiNW with many 
randomly positioned defects puts strong requirements on 
the method. The quasi onc-dimcnsional nature of the 
SiNWs requires on the one hand an atomistic model tak- 
ing quantum effects and charge transfer around the de- 
fect into account. On the other hand the method should 
be able to treat more than 10^ atoms and include many 
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scattering events due to the /zm length of the wires. 

In this work, two methods are used to study the eflFect 
of disorder including many randomly placed H- vacancies 
in H-passivated long SiNWs. Both methods are based on 
ab initio calculations and scale linearly with the length of 
the sample. The first approach uses a relatively new real- 
space method, developed by S. Roche and D. Mayou and 
co-workers over the last decade to study transport prop- 
erties primarily in carbon nanotubes and quasi-periodic 
systemsi^^i^°i^^i^^i^^i^'^i^^i^^i^'^ The method is based on 
the Kubo-Greenwood formalism22iSS rewritten in a real- 
space framework, and we will refer to it as the Kubo 
method. The second and more well-known approach is 
based on the Landauer formula where the conductance 
is found by recursive calculations of Green's functions 
(GFs). We will refer to this as the GF method. 

The Kubo method was shown to predict the elastic 
MFP at the Fermi energy in randomly disordered car- 
bon nanotubes (CNTs)i^, in agreement with a Fermi's 
golden rule estimate^. Besides that, we are aware of no 
comparison with other theoretical methods. One of the 
primary goals of this paper is to report such a compar- 
ison over for several energies. We show that the Kubo 
and GF methods are in general in qualitative agreement, 
however, with quantitative differences especially around 
band edges. We analyze the pros and cons of the two 
methods and give an assessment of when they should be 
applied and which quantities they can calculate. 

The rest of the paper is organized as follows. In sec- 
tion^jwe summarize the two numerical methods and de- 
scribe how a Hamiltonian for a long SiNW is constructed 
from ab initio calculations. Results concerning both un- 
passivated as well as hydrogen passivated SiNWs are pre- 
sented in section Unl We end up with a discussion of the 
applied methods and the results in section Hvl 



II. METHODS 
A. Real space Kubo method 

In the real space Kubo method, which is derived from 
the Kubo-Greenwood formula, transport properties are 
determined by calculating the time propagation of wave 
packets in real space. The central quantity is the time- 
and energy dependent diffusion coefficient, V{E,t), de- 
fined by 



l Tr{(X(t)-X(0))^,5(i;-H)} 
^(^''^ - 1 Tr{<5(£;-H)} ' 



(1) 



where X(t) is the position operator along the wire direc- 
tion written in the Heisenberg representation, H is the 
Hamiltonian matrix and the trace Tr{(5(_E — H)} is the 
total electronic density of states (DOS). Following Trio- 
zon et alm^^ an efficient evaluation of the traces can be 
carried out by using a relative modest number (< 10) of 
random phase states, {ipr)- The coefficients of Itjjr) which 



are labeled from —N/2 to N/2, with TV being the total 
number of orbitals in the wire, are initially non-zero only 
in the central part of the nanowire, i.e. 




o2iTra(j) 



for - Nr/2 <j< Nr/2, 
otherwise. 



(2) 

where a{j) is a random number in the interval [0, 1[ and 
typically Nr ~ 10^. The number of random phase states 
needed to accurately estimate the traces is not known a 
priori and the convergence of the results must be checked. 

Time evolution of the random phase states can be ef- 
ficiently computed by expanding the time evolution op- 
erator, e~*^*/'^ in the orthogonal set of Chcbyshcv poly- 
nomials. Each term in the traces in Eq. is a local 
density of state which is calculated using a continue frac- 
tion technique. This is the most time consuming part 
of the calculations and involves ^ 10^ operations with 
the Hamiltonian for the considered systems to resolve 
the closely lying energy bands. The convergence of the 
continued fraction scheme must be separately verified. 

The elastic MFP, le, is found from^^. 



le{E) = 



ma.x{V{E,t), t > 0} 



(3) 



where v{E) is the group velocity at energy E. 

Notice that there is no requirement that the Hamilto- 
nian is periodic. Nor are there any leads that connect to 
a device area. We also stress that the time it takes to 
calculate the diffusion coefficient at many energies is not 
much longer than for a single energy. The reason for this 
is that the random phase states contain all energy com- 
ponents and the time evolution is energy independent. 
Moreover, the primary numerical task in the continued 
fraction scheme is a mapping of the original Hamiltonian 
to a smaller tridiagonal matrix. This mapping is energy 
independent, and it is therefore relatively fast to compute 
the MFP for the whole energy spectrum. 



B. Recursive Green's function method 

The second numerical method we have applied is based 
on the Landauer formalism described in detail in e.g. Ref. 
Is^ The general setup is illustrated in Fig. A device 
region (D) is connected to a left (L) and right (R) semi- 
infinite lead. The device area consists of M sub-cells and 
is described by the Hamiltonian H^i 
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(4) 



The sub-cells are chosen so large that only nearest neigh- 
bour cells couple. Generally, the sub-cells do not need to 
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have the same size. The leads described by and Hn 
respectively are assumed to have a semi-infinite structure 
consisting of equal unit cells, Hq. The coupling matrices 
between the leads and the device area are denoted 
and Vij. The Hamiltonians are calculated with a non- 
orthogonal basis set (see section III (1 such that for each 
Hamiltonian matrix we also have a corresponding overlap 
matrix S. 




Vl 



D 



Vr 



R 



FIG. 1: The device region (D) is divided into M sub-cells. 
The sub-cells are so large that they only couple to the nearest 
neighbour cells. Also, the left (L) and right (R) leads consist- 
ing of equal unit cells, Ho, couple only to the first and last 
cell of the device region respectively. 



To calculate the length dependent conductance of the 
wire, we initially find the surface Green's functions, 
Glj^{E) = (-BSl,^ - HL,fl)~\ of the isolated leads by 
a standard decimation procedurei^^ The device area is 
subsequently grown by adding one sub-cell at a time and 
calculating the Green's function 



(■0 



(5) 



where Hii{E) describes the coupling to right lead and 
is defined through G"j^{E) = {ESq - Hq - SJ^(^;))~^ 
S L (E) also takes the coupling to the rest of the device 
area into account and is calculated from the previous 
Green's function as = V^Li,,Gi_i(£;)V,_i,,. 

In each growth step we calculate the conductance of 
the wire consisting of the first i sub-cells as 



g{E,U) = ^Tr [gI{E) Tr{E) G,{E)Tl{E) 



(6) 



where Li is the length of the grown device region, 
Tji,l{E) = — 2 Im [Sfl.L(£')] and where the trace is per- 
formed over the states in the device region. 

Sample averaging is performed over 200 different con- 
figurations giving a mean conductance, (g). The cor- 
responding resistance is found as R = 1/(5)- For wire 
length L <ti ^, with ^ being the localization length, the 
resistance increases linearly as R{L) = Rq -\- RoL/le, 
defining the MFP, In the localization regime the 

resistance increases exponentially and we calculate the 
localization length as^ 



lim 



2L 



L^oc (In g) ' 

where the limit is taken such that ^ is converged. 
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C. Constructing the Hamiltonian from ab initio 
calculations 

The atomic and electronic structure of the SiNWs is 
found from ab initio calculations using the density func- 
tional theory (DFT) package SiestAi^ Using first prin- 
ciple calculations it is relatively straightforward to in- 
troduce various defects such as vacancies, adatoms or 
dopants. This is not the case when using standard tight- 
binding parameters. 

We have used a minimal single-^ basis set'''', with 
4 orbitals (one 3s and three 3p) on the Si atoms and 
one on the H atoms, to represent the one-electron wave 
function. The minimal basis set is applied in order to 
speed up the subsequent transport calculations. We have 
used norm-conserving pseudopotentials of the TrouUier- 
Martins typeS and the Generalized Gradient Approx- 
imationS^ for the exchange-correlation functional. The 
calculations are performed on super-cells containing five 
wire unit cells (see Fig. The reciprocal space has 

been sampled with a converged grid of I x I x 2 fc points 
following the Monkhorst-Pack scheme.^^ 

When modelling a defect this is placed in the middle 
of the supercell to ensure that the effect of the defect 
is confined to the central region and is not affected by 
the periodicity and the inter-cell coupling terms Vi^i in 
Eq.(4) are independent on the specific type of defect. 
The atomic postions of all the atoms in the cell contain- 
ing the defect and in the two neighboring cells have been 
fully relaxed, until the maximum force was smaller than 
0.04 eV/A. This is an important point, because the lo- 
cal distortion induced by the defect can have important 
consequences on its scattering properties. 

To model a long (L > 100 nm) SiNW with a random 
distribution of defects (section If f f fj|l we perform a full 
Siesta calculation for a pristine wire and for each type 
of defect. The wire is grown by adding pieces from the 
different calculations ending up with a wire structure as 
in Fig. n]and a Hamiltonian as Eq. (gj. The sub-cells 

H^-* can either be pristine parts or contain a defect. 

Usually when joining pieces from different calculations 
the Fermi energies should be aligned to ensure charge 
conservation. However, for the SiNWs with hydrogen 
vacancies, a dangling bond (DB) state forms in the band 
gap and pins the Fermi energy causing a large shift com- 
pared to a pristine wire. Alignment of the Fermi energies 
would therefore lead to an unphysical shift in energies 
also for sites far away from the vacancy. We therefore 
align the Si 3s on-site energy instead, taking as reference 
the first Si atom in the supercell. 

The minimal basis set implies that we can not expect 
the conduction bands to be accurately described, and we 
will thus only focus on the valence bands. Furthermore, 
the Hamiltonian is truncated by removing the smallest el- 
ements, while keeping the band structure constant with 
a tolerance of 10~^ eV. This leads to a reduction of non- 
zero matrix elements of around 80%, and we thus effec- 
tively end up with a tight-binding-like Hamiltonian. 
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Siesta uses a non-orthogonal atomic basis set giving 
rise to an overlap matrix S. The Kubo method, how- 
ever, requires an orthogonal basis set and the Siesta 
output cannot therefore be used directly for which reason 
a Lowdin transformation is performed. For each Siesta 
calculation the Hamiltonian, H, containing five unit cells 
described in the non-orthogonal atomic basis, is mapped 
to a new one, H, with an orthogonal basis: 

H ^ H = S-i/2fjs-i/^ (8) 

In the same way as described above, a long wire is built 
by extracting parts from both pristine and defected or- 
thogonalized Hamiltonians and joining the pieces. The 
Lowdin transformation leads to a longer range of the ba- 
sis orbitals and thus more non-zero elements in the trun- 
cated Hamiltonian. This in turn implies longer calcula- 
tion times. 



D. Fermi's golden rule 

We next compare the two numerical approaches with 
results obtained using Fermi's golden rule (FGR). We 
consider scattering due to the scattering potential V be- 
tween unperturbed Bloch states, \n,k), of the pristine 
wire. The transport relaxation rate, 1/t„, from band n 
at energy E is found as 

J.,,, = ^ / dk'J2\{k',m\V\n,kW 

Tn{En{k)) h J ^ 

x{l - COS Okk') S{E„,{k') ~ En{k)) 
= -^Y.\^-k"MV\n,kWn^{E),{%) 

m 

where the summation is over bands, n^{E) is the DOS 
corresponding to band m and the energy of the final state 
fulfills Em{—k") = En{k) = E. 9kk' is the angle between 
the initial and final wavevector. Only scattering from a 
forward to a backward propagating state contributes to 
the rate, i.e. dkk' = ^tt yielding a factor 2. The per- 
turbing potential, V , is for the Anderson on-site disorder 
f section IIII A|) simply a diagonal matrix with the diago- 
nal elements being random numbers in a given interval 
(see below), and the bar denotes an average over different 
configurations of V . The MFP for electrons in band n is 
calculated as 

ln{E)^Vr,{E)T,,{E), (10) 

where Vn{E) is the group velocity at energy E in the 
n'th band. We calculate the total MFP as the mean 
value from the individual contributions from the bands: 

^e(i?) = ^E'"(^)' (11) 

where N{E) is the number of bands at energy E. 



III. RESULTS 

A. Surface reconstructed wires 

The first structures we consider are un-passivated, sur- 
face reconstructed SiNWs as illustrated in Fig. |2] show- 
ing a side-view (top) and a cross-section view (bottom 
left) of a wire containing 5 unit cells, separated by the 
red dashed lines in the upper panel. This structure was 
recently studied by Rurali et alw^ using DFT calcula- 
tions. There are 57 atoms in one unit cell, the diame- 
ter is ^1.5 nm, and the length of the unit cell in the 
growth direction is 0.55 nm. Notice that the structure 
at the surface differs significantly from the bulk of the 
wire. Ref. 12 showed that this particular surface struc- 
ture makes the wire metallic, which is evident in Fig. |21 
(bottom right) showing the band structure around the 
Fermi energy (marked by the dashed line). Four bands 
are crossing the Fermi level (two being degenerate). 




FIG. 2: (Color online). Top: Side-view of the un-passivated 
SiNW containing 5 unit cells, separated by the red dashed 
lines. Bottom left: End view of the un-passivated SiNW. Bot- 
tom right: Band structure around the Fermi energy (dashed 
line). Four bands are crossing the Fermi energy, two being 
degenerate (with E = — O.OSeV at the F— point). 

To further investigate the metallic surface states and 
to compare the Kubo and the GF method we add random 
on-site noise to either the surface atoms only or to the 
bulk atoms only. The relevant diagonal elements in the 
Hamiltonian are thus changed according to 

H,, -^H,,+,5,;, (12) 

where 5i takes values with equal probability in the inter- 
val [— Ae/2; Ae/2] with Ae being the disorder strength. 
This Anderson model for electronic disorder is simple 
and widely applied, very recently Zhong et al^ used 
it to model surface disorder in shell-doped nanowires. 
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Whether this simple model adequately describes phys- 
ical defects merits a separate discussion, given at the end 
of section IIII Bl 

Fig. O shows T>{Ef, t) revealing that the bulk disorder 
has a very small effect on the transport properties yield- 
ing an almost linearly increasing diffusion coefficient, a 
characteristic feature of ballistic transport. The surface 
disorder, on the other hand, leads to diffusive transport 
with V{E,t) almost constant for t > 150 fs. The pro- 
nounced differences of the surface- and bulk-disordered 
systems show that the conduction around the Fermi en- 
ergy almost entirely takes place along the surface atoms, 
in qualitative agreement with the conclusions drawn in 
Ref. [13 
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FIG. 3: (Color online). Time dependent diffusion coefficient 
at -E = Ef- The bulk disorder has little effect and the trans- 
port is ballistic. The edge disordered system shows diffusive 
behaviour with a constant diffusion coefficient for t > 150 fs. 

Fig. ^ (upper panel) shows the MFP, le{E), vs. energy 
calculated using both the Kubo method (soHd blue) and 
the recursive GF approach (black squares). The dashed 
red line is an analytical estimate obtained using FGR. 
The disorder strength is Ae = 0.4 eV and only the sur- 
face atoms are perturbed. The Kubo and FGR results 
are obtained with a fine energy resolution whereas the 
GF results are calculated only for relatively few, discrete 
energies, illustrating the advantages in the Kubo method 
in calculating properties at many energies in one calcu- 
lation. 

Generally all three methods agree qualitatively and in 
the interval —0.1 ^ E < 0.2 the results are even quan- 
titatively consistent. In this energy range several bands 
exist, see Fig. |21 (bottom right) and thus there are more 
possible back scattering processes giving a larger scatter- 
ing rate and thus a shorter MFP. Notice, that lowering 
the mean free path does not necessarily mean reducing 
the conductance, because we also have more charge car- 
rying states. 

Around the band edges at E ^ -0.09 eV the MFP 
calculated with FGR drops sharply whereas the values 
obtained with the GF and Kubo methods drop more 
slowly. This difference is caused by the relatively large 
disorder strength, Ae = 0.4 eV, which broadens the DOS 
and smears out the sharp features. For smaller disorder 
strengths the GF results resemble the FGR values more. 




-0.1 0.1 
Energy (eV) 



E = -0.2 eV 
□ E = 0.0 eV 
V E = 0.25eV 




FIG. 4: (Color online). Top: Mean free path, le, vs. en- 
ergy for Ae = 0.4 eV. The solid lines are obtained with the 
Kubo method, the GF results are marked by squares while 
the dashed lines shows the FGR results. GF results are av- 
erage values for 200 different samples while the Kubo results 
are mean values of 10 different samples. Bottom: Scaling of le 
vs. 1/Ae shown on a logarithmic scale. Circles, squares and 
triangles are calculated with GF while the lines are obtained 
using FGR at the same energies. 



which is illustrated in the lower panel in Fig. 0] The 
figure shows, on a logarithmic scale, the MFP vs. in- 
verse disorder strength, 1/Ae, at three different energies. 
The points are obtained with the GF method, while the 
lines are obtained using FGR. For weak disorder, the GF 
results scales as le oc (1/Ae)'^ in accordance with FGR, 
whereas the GF results for strong disordered systems de- 
viate from the (1/Ae)^ dependence. In this regime the 
first order perturbation applied in FGR does not fully 
suffices and reliable results must be obtained with more 
elaborate approaches such as the GF method. 

Around 0.24 eV, the Kubo method fails to find the pro- 
nounced peak obtained with both the GF method and 
with FGR. This difference is probably caused by two ef- 
fects. The first reason is again a broadened DOS, since 
the GF results show a similar deviation from the FGR 
values as seen around E = — 0.09eV. The second and 
more important reason to the differences is due to the 
numerical inaccuracies in calculating the density of states 
in the Kubo method. The inaccuracies are especially im- 
portant around van Hove singularities at the sub-band 
edges. E ~ 0.24 eV marks the band edge for the two 
degenerate bands, and due to the finiteness of the sys- 
tem considered in the Kubo calculation, the van Hove 
peaks in the DOS will unavoidably have decaying tails 
at larger energies. This implies that the calculated den- 
sity of states will be too large causing ^{E, t) and thus 
le to be correspondingly smaller. The energy separation 
between the two sub-band edges around E = 0.24 eV is 



6 



only ~ 0.05 eV, which is less than 0.3% of the total band- 
width, W 20 eV. This makes it numerically difficult to 
resolve the detailed features with the continued fraction 
technique used by the Kubo method. 



B. Passivated wires 

The surface reconstructed SiNWs are both physically 
and technologically exciting, but probably also very frag- 
ile objects, since small changes in the surface such as de- 
fects or adatoms presumably will change the performance 
drastically. Moreover, the wires produced are most often 
surface passivated by either SiO^; or hydrogen, and the 
focus in this section will be on such wires. The passivated 
wires are semiconducting, often with a direct band gap 
that increases for small diametersSiiSii^. The simplest 
way to model surface passivated wires is to add hydro- 
gen atoms to the surface such that all the Si dangling 
bonds are passivated. Such wires resemble qualitatively 
those reported by Ma et al.^. 

Figure [SI shows the cross section of the wire (left). The 
unit cell contains 57 Si atoms and 36 H atoms labeled 
with a number from 1 to 36 as indicated in Fig. [Sl(left). 
The band gap is found to be 2.84 eV4i 




FIG. 5: (Color online). Left: Cross section of the H passivated 
wire. Right; Energy dependent conductance of a pristine wire 
and of wires containing a single vacancy of number 1-3. 

To investigate the influence of surface defects we in- 
troduce hydrogen vacancies. The vacancies are labeled 
correspondingly to the removed hydrogen atoms. Note, 
that due to symmetry there are only five topologically 
different vacancies. The conductances of wires with only 
a single vacancy, corresponding to one of the H atoms 
1-3 being removed, are shown in Fig. |S1 (right) for ener- 
gies close to the valence band edge {E — OeV). Clearly, 
vacancy 1 and 2 scatter the most, the reason probably 
being that the wave function (for the pristine wire) has a 
pronounced larger weight on these H atoms. Notice, for 
energies —0.15 < E < — 0.05eV, one channel is almost 
completely closed by vacancy 1. The vacancies 3-5 give 
almost the same conductances. 

We model a wire with randomly missing hydrogen 
atoms by performing a Siesta calculation for each possi- 
ble vacancy position (one of the 36 H atoms is removed) 
and adding pieces from the different calculations. We 
measure the vacancy concentration by the average dis- 



tance, {dn), in the wire direction between two vacancies. 
Each unit cell can only contain one vacancy thus setting a 
lower limit for (du) at the unit cell length, a — 0.56 nm. 
The MFP for (dn) = 5.6 nm calculated with the GF 
method is shown in Fig. |S1 (solid) revealing a strong en- 
ergy dependence. In the interval —0.15 < E < — 0.05eV, 
we find le ~ 50 nm, while for energies around — 0.35eV 
the MFP in on the order of 1 ^m. Comparing with esti- 
mated phonon scattering MFPs of more than 500 nm^i^ 
the application of the elastic scattering model applied in 
this work is justified for most of the energies. Moreover, 
at several energies the vacancy scattering might be the 
dominant even at room temperature. However, around 
the peak at E = — 0.4eV where the calculated MFPs 
exceeds 1 /im, other scattering sources are likely to dom- 
inate. 

Assuming that all bands at a given energy have 
the same reflection probability, Ri{E) = {To{E) — 
Ti{E))/To{E), where Tq is the total transmission of a 
pristine wire, and Ti is the total transmission through a 
wire containing a single vacancy with number i (shown in 
Fig. El (left)), the MFP in a wire with only vacancies of 

type i, can be estimated as li^\E) = {dn) / Ri{E). The 
total MFP can be estimated using Matthiessen's rule, 
such that f^^ — X]i(^e*'')~^, and the result is shown in 
Fig. IHlas the dash-dotted line. It is evident that the sim- 
ple estimate, which ignores interference effects between 
successive scatterers, accurately reproduces the results 
found by sample averaging over vacancy configurations. 

The Kubo method requires a Hamiltonian describing 
a wire that is longer than the largest mean free path in 
order to avoid boundary effects. As seen from Fig. 
we therefore need a wire of length L > 1500 nm consist- 
ing of more than 3,000 unit cells and thus N > 8 ■ 10^ 
orbitals. In our current implementation this causes com- 
puter memory problems, and we have not been able to 
obtain reliable results with the Kubo method for the va- 
cancy scattering. 
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FIG. 6: (Color online). Mean free path vs. energy for a con- 
centration of vacancies corresponding to (dn) = 5.6 nm. Solid 
black line corresponds to wires containing all possible vacan- 
cies. Dashed-dotted blue is obtained using Matthiessen's rule 
for single scattering events and the dashed red curve corre- 
sponds to a wire with Anderson disorder (Ae — 1.3 eV). 

We next examine whether the vacancies effectively can 
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be modelled by adding Anderson disorder. The calcu- 
lated MFP for a surface disordered wire, i.e. with no 
vacancies but rather on-site disorder at all orbitals at 
the surface Si atoms, is shown in Fig. El (dashed red). 

For a disorder strength Ae = 1.3 eV the Anderson 
model fits the vacancy results at energies E > — 0.35eV, 
although the peak around E — —0.05 eV is much less pro- 
nounced. The small shift of the peaks around E = —0.23 
eV and E = —0.05 eV is due to a broadened DOS in the 
Anderson disordered wires. For energies below —0.35 eV, 
the Anderson model deviates significantly from the va- 
cancy results. The pronounced peak in MFP is not found 
in the Anderson model which gives almost constant val- 
ues up to a factor of three lower than the vacancy results. 
Besides a broadened DOS the differences probably arise 
because the Anderson model is a too simple model not 
capturing all the physics. Since the Anderson disorder 
only reproduces the vacancy results at some energies we 
conclude that the effects of vacancies cannot accurately 
be modelled by the simple on-site disorder. Moreover, 
the value of the disorder strength, Ae, has no clear con- 
nection to an actual vacancy concentration. 



lower than e.g. the MFP of le 
energies around E = — 0.35eV. 



1 /im at slightly lower 




50 100 150 200 250 300 
Length (nm) 

FIG. 7: (Color online). Length dependent resistance for at the 
energies E = — 0.3eV (dashed blue), E = — 0.15eV (dashed- 
dotted green) and E = — O.OSeV (solid red). The average 
distance between vacancies is {(Ih) = 2.8 nm. The inset shows 
the scaling of le vs. average distance, {(Ih), between vacancies 
at the same three energies as above. 

Fig. [7|shows the length dependent resistance for three 
different energies at vacancy concentration corresponding 
to (dn) ~ 2.8 nm. For length L < 200 nm the resistance 
increases linearly at all energies, with the slope determin- 
ing the MFP. The inset shows the scaling of le vs. {(Ih) 
at the same three energies as in the main figure. The 
points are calculated with the GF method and the lines 
are linear fits, clearly revealing a linear relationship be- 
tween the MFP and the average inter-vacancy distance. 
For length L > 200 nm the resistance corresponding to 
E = — 0.15eV (dash-dotted green) starts to increase ex- 
ponentially thus entering the localization regime. The 
localization length at this energy is ^ = 110 nm, much 



IV. DISCUSSION 

In this paper we have studied electronic transport in 
SiNWs and calculated the influence of disorder on the 
mean free path. Our model is subject to a number 
of limitations and approximations. We apply a single- 
electron model and consider the linear response regime. 
Moreover, the minimal basis set may limit the accuracy. 
Also the spin orbit coupling is not included which is nec- 
essary to describe the details around the valence band 
edge in bulk Si. Compared to the experimentally realized 
SiNWs, the wires considered here are quite thin although 
comparable to the wires reported in Ref. U The SiNW 
structures used in this paper have rounded, rather than 
perfectly square angles. It has been proposed in the liter- 



aturi 



.44.45.46 



that smooth angles would be favored during 
the growth process with respect to the sharp angles that 
naturally arise from the square symmetry of the [100] 
cleavage plane. The topic has been discussed at some 
details elsewhere^^ for surface reconstructed wires, con- 
cluding that at nanometric diameters no clear difference 
emerges. At the same time the electronic structure seems 
to be only marginally affected. Since all calculations in 
this work are performed on relaxed structures fully based 
on ab initio calculations without any use of fitting param- 
eters, we expect, in spite of all the limitations, to capture 
the correct trends in the transport characteristics. 



A. Methods 

Two numerical methods were applied and compared 
to each other: A real-space Kubo method, and a re- 
cursive Green's function method. The two approaches 
each have their advantages: In calculating the MFP at 
many energies, the Kubo method is advantageous, since 
the diffusion is readily found at many energies in a single 
calculation, whereas the GF method requires a full calcu- 
lation at each energy. If the focus is on a few energies but 
many different disorders, the GF method is the preferred 
choice. For metallic systems, where one mainly is inter- 
ested in the properties around the Fermi level, the GF 
method thus seems to be the method of choice - the par- 
allel computation of many energies in the Kubo method 
is not needed. The Kubo method seems more applicable 
to semiconducting systems, since it is physically possi- 
ble to move the Fermi level with a gate voltage, thereby 
scanning several energies. 

The Kubo method requires a Hamiltonian describing a 
wire that is longer than the largest MFP. For the weakly 
disordered wires with long MFPs we need wires of length 
L > 1 /im containing more than 10^ atoms resulting in 
very large matrices. In our current implementation this 
causes memory problems and the Kubo method failed to 
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converge for the H-passivated SiNWs. The GF method, 
on the other hand, involves only calculations with the 
small sub-cell Hamiltonians and it suffices to consider 
wires grown to a length L « 50 nm to get an accurate 
estimate of the initially linear resistance vs. length curve 
(cf. Fig. 13). Moreover, it proved to be numerically dif- 
ficult to resolve the detailed features in the energy spec- 
trum with the Kubo method, which led to erroneous re- 
sults near band edges. The difficulties arise because the 
Kubo method does not take any semi-infinite periodic 
leads into account as in the GF method. There, the DOS 
is readily calculated to arbitrary accuracy from the sur- 
face Green's function by using the periodic structure 
of the leads. However, the absence of periodic leads in 
the Kubo method can also be a great advantage since it 
allows to study non-periodic systems such as incommen- 
surable multi-walled carbon nanotubesii^ 

The growth procedure in the GF method involves in- 
versions of the relatively small matrices describing the 
sub-cells. For thin wires as those considered in this work, 
with N ~ 250 orbitals in each unit cell, the inversion 
step is not a critical issue. However, for thicker wires the 
number of orbitals increases quadratically, and due to an 
0{N^) scaling of the inversion step, the whole procedure 
scales as 0{d^), with d being the wire diameter. Con- 
trary, the Kubo method scales linearly with the number 
of orbitals, and thus as 0{d^). 

For systems with a periodic structure, as the SiNWs 
considered here, we generally find the GF method to be 
the preferred choice, given that a rough energy resolution 
is sufficient. Thicker wires with relatively short MFFs 
favor the use of the Kubo approach. 



adding the individually calculated MFPs to an effective 
MFP by applying Matthiessen's rule. However, it is not 
obvious that a similar rule applies for the localization 
lengths thus emphasizing the importance of considering 
systems containing many differently placed vacancies as 
opposed to only considering single scattering events. It 
was further shown that an attempt to model the vacan- 
cies with an effective Anderson disorder gave satisfactory 
results only in a limited energy interval. We therefore 
conclude that the Anderson model is a too simplistic de- 
scription of the real vacancies. 

The MFP was shown to be strongly energy depen- 
dent, and for relatively long wires the resistance can be 
changed by orders of magnitude within a 0.1 eV shift of 
the Fermi energy thus causing a transition from the dif- 
fusive (Ohmic) regime to the localization regime. The 
strong energy dependence might be utilized in sensor ap- 
plications where the presence of a single virus acts as a 
local gate shifting the Fermi energy^ This could possi- 
bly cause a transition from the Ohmic to the localization 
regime thus changing the resistance of the wire dramati- 
cally. However, more careful work has to be done before 
firm conclusions can be stated. 

For relatively strong disordered wires, the MFP is well 
below 500 nm for a large energy range. Compared to es- 
timated phonon scattering MFPs of more than 500 nmi^, 
the elastic scattering model applied in this work is justi- 
fied. Comparing the results obtained in this work with 
the estimated long phonon MFP we suggest that impu- 
rity and defect scattering could be the dominant scatter- 
ing source even at room temperature. 



B. Results 

In un-passivated, surface reconstructed wires, Ander- 
son disorder was added to the surface atoms affecting the 
transport properties significantly. Disorder in the bulk 
had, on the other hand, no significant influence. 

In hydrogen passivated wires surface disorder was in- 
troduced by randomly removed hydrogen atoms. We find 
that it suffices to consider single scattering events and 
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